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ABSTRACT 

In  the  area  of  ranking  and  selection,  recent  papers  have  presented 

procedures  which  attempt  to  allow  experimenters  to  take  at  least  partial 

advantage  of  more-favorable  configurations  of  the  population  parameters 

(compared  to  the  least -favorable  configuration).  The  current  paper  presents 

a  computer  program  which  allows  experimenters  to  dispense  with  all  of 
inequalities,  an  indifference -zone,  and  known  equal  variabilities  in  their 
assessment  of  the  probability  of  correct  selection.  It  has  additional  uses 
in  implementation  of  recently -proposed  sequential  adaptive  selection  procedures 
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ESTIMATED!!  OF  THE  P(CS)  WITH  A  COMPUTER  PROGRAM 


by 

Edward  J.  Dudewicz 


1.  INTRODUCTION 

When  an  experimenter  uses  the  indifference-zone  ranking-and-selection 
procedures  oioneered  by  Bechhofer  (195*0  and  recently  explicated  by  Gibbons, 

Olkin,  and  Sobel  (I977)j  be  sets  his  sample  size(s)  to  guarantee  a  P(CS) 
(probability  of  correct  selection)  at  least  P*  whenever  a  minimal  separation 
(of  the  best  population  from  the  other  populations)  of  nt  least  o*  exists; 

a  separation  of  6*  is  then  least -favorable.  The  actual  separation  is  likely 

to  be  more  favorable  than  least -favorable  (and  the  P(CS)  thus  >  P*) ,  and 

experimenters  would  like  to  take  advantage  of  this  favorableness  after  the 

experiment  has  been  run.  Previous  approaches  to  this  problem  by  Olkin, 

Sobel,  and  Tong  (1976)  and  Anderson,  Bishop,  and  Dudewicz  (1977)  all  involved 
inequalities.  In  this  paper  we  give  a  computer  program  which  allows  assess¬ 
ment  of  the  P(CS)  without  use  of  inequalities,  is  simple  and  easily  implemented 

on  even  small  computers,  and  does  not  use  an  indifference-zone.  It  is  also 

directly  useful  when  (as  is  often  the  case)  the  samples  are  drawn  before  the 

statistician  is  consulted,  and  he  must  evaluate  their  adequacy.  It  has 

additional  uses  in  implementation  of  recently  proposed  sequential  adaptive 

selection  procedures  of  Tong  (1978). 
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While  our  approach  is  oriented  towards  the  practitioner  and  computer 
software  provider, -in  this  paper,  it  should  he  noted  that  the  area  is  also 
one  of  recent  intensive  theoretical  investigation  by  Olkin,  Sobel,  and  Tong 
(1979)>  Bofinger  ( 1900) ,  and  .  *  cm  (1980  a,  b) . 
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Note  that  Olkin,  Sobel,  and  Tong  (1979)  consider  only  the  case  of  o 
known.  Work  of  Bofinger  (1980)  contradicts  their  result  on  asymptotic 
normality  of  their  ?(cs)  estimator.  Falti.i  (1980a)  shows  this  estimator 
is  biased  on  the  high  side  when  k  ■  2,  while  Faltin  ( 1980b)  gives  an 
alternate  quantile  unbiased  estimator  (also  for  the  case  k  =  2) .  By 
contrast,  our  estimator  allows  for  unknown  and  heteroscedastic  variances  and 
does  not  require  either  bounds  or  tables:  the  experimenter  makes  one 
inexpensive  run  of  a  simple  computer  program  (presented  in  full  detail  with 
numerical  examples  below) . 
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2.  TVE  P(CS)  FUNCTION 

If  the  experimenter  has  obtained  a  random  sample  of  size  from 

population  which  yields  normally  distributed  observations  with 

a 

unknown  mean  and  variance  ^  ,  1  <  i  <  k,  and  wishes  to  evaluate  the 
adequacy  of  these  sample  sizes  when  the  means  procedure  ("select  the 
population  yielding  the  largest  sample  mean  a3  having  the  largest 

population  mean")  is  used  (and  all  n^+ . . .+  n^  observations  are  independent:* 


he  will  need  to  assess  the  function 
®  k-1  _  , 


p(cs)=  ;  {  *  z  + 

1  i-1  v  a(i)  v  n(k) 


M[k]  ~  ^[13  )} 


)}  4  (z)dz 


where:  <  ...  <  uj-^j  denote  uk  in  numerical  order;  a^,  n^, 

are  the  standard  deviation  and  sample  size  of  the  population  with  true  mean 
4j- 1] ;  t,  4>  are  the  standard  normal  distribution  function  and  density  function. 

Note  that,  with  error  5  (0  <  S  <  .005  since  *  (2.8l)  =  .99T5)> 

2.81  k-1  . - 

P(CS)  -  f  (  IT  «  (-isi  —til  Z  +  )}  f(z)&z. 

Once  the  pj^j,  are  estimated  from  the  samples,  this  can  be  estimated 


simply  by  Monte  Carlo  methods  on  even  small  computers,  with  ease  compared  to 

quadrature  methods,  and  with  accuracy  compared  to  methods  involving  inequalities. 
A  computer  program  is  given  in  Section  3,  with  an  example  in  Section  h. 
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COMPUTER  PROGRAM 


The  following  computer  program  reads  K  (the  number  of  populations'  and 

S(t),  II (l)  (the  sample  mean,  sample  standard  deviation,  and  sample 

size  of  population  T),  1  <  I  <  K,  from  data  cards  and  then  uses  '.0,000 

)'ionte  Carlo  trials  to  estimate  ?(CS  with  a  standard  error  not  exceeding 

.005/2,  so  that  (with  Jo*  confidence,  also  taking  account  of  the  truncation 
error  €)  we  estimate  the  Pecs'!  integral  within  .01.  A  numerical  example 
io  given  in  section  4. 
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non  n  oonn  nnnnno  anftrtftonnnoo 


f 

THE  FOLLOWING  COMPUTER  PROGRAM! 

1.  READS  K  (THE  NUMBER  OF  POPULATIONS)  FROM  A  DATA  CARD; 

2.  READS  XM(  I) ,S< I) ,N< I)  (TOE  SAMPLE  MEAN,  STANDARD  DEVIATION 

AND  SIZE  OF  POPULATION  I),  1<=I<=K,  FROM  K  DATA  CARDS; 

3.  USES  40,000  MONTE  CARLO  TRIALS  TO  ESTIMATE  THE  P(CS)  WITH 

A  STANDARD  ERROR  NOT  EXCEEDING  .003/2,  SO  THAT  (WITH  93 
PERCENT  CONFIDENCE,  ALSO  TAKING  ACCOUNT  OF  THE 
THTJUCATION  ERROR)  WE  ESTIMATE  THE  P(CS)  INTEGRAL  WITHIN 
.01. 

INTEGER  R, SEED 

DIMENSION  AiLZQ)  ,S(50)  ,N(GO)  ,R(50) 

PERSONS  WITH  K>50  SHOULD  INCREASE  ALL  FOUR  DIMENSIONS  ABOVE 
TO  THAT  NUMBER  BEFORE  RUNNING  THIS  PROCRAM.  THEY  SHOULD 
ALSO  INCREASE  THE  DIMENSIONS  OF  YM(  . )  AND  S( . )  IN  SUB¬ 
ROUTINE  VSORT. 

READ! 3. 20)  K 

20  FORMAT! 110) 

DO  30  I  =  1 ,  K 

READ! 3,21)  XM(  I) ,S(  I) ,N( I) 

21  FORMAT!  F 10. 0,F 10.0. 110) 

30  CONTINUE 

WE  NOW  FIND  INTEGERS  R(  1) . R(K)  SUCH  THAT 

XM<  ll(  1 ) )  <  =XM(  R(  2) )  <  =  .  .  .  <  =XM(  R(  K)  )  BY  A  VIRTUAL  SHELL  SORT. 

CALL  VSORTf  XM, R, K) 


SELD= 067021 
1  =  0 

JCNT=0 


WE  NOW  ENTER  THE  MONTE  CARLO  LOOP. 

31  CONTINUE 

1F( I . CE. 40000) GO  TO  34 
GALL  URN20(  SEED ,  RANX,  1 ) 

CALL  IJRN20'  SEED ,  RANY.  1 ) 

RAHX=  ( 11ANX-0 . 3 )  *2 . 0*2 . 0 1 

11AHY=  RANY*(  1 . 0/SQRT(  2 .  ©*3 .  14159)  ) 

PROD*  1 . 0 

■<K  =  K-l 

1)0  32  L=1,TCK 

TERM*  XM<  IK  !C) )  -XM(  R<  L) ) 

TERM®  TERM/S ( IK  L.) ) 

AN  =  FLOAT!  H  .iKL))) 

TERM  =  TERMS  J'AilT!  AN) 

FACT=S(  R(  K) )  /S(  IU  L) ) 

CN  =  FLOAT!  N!  IK  L) ) ) 

ON  =  FLOAT!  N<  K<  K)  )  ) 

BN  =  CN/DN 

FACT  =  FACY*'JLT.<  ON; 

FACT*  FACT*  RAN  X+  TERM 
TERM*  DCDFNt  FACT) 
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finrififi 


PROD  =  PROD*TCRM 

32  CONTINUE 

FN=PR0D/'3QJIT!2.0*3. 14159) 

FM=> FN*EXF(  -0 . 3*RANX*RANX) 

IF! FN.LE. RAN?) CO  TO  33 
JCNT= JCNT+ 1 

33  CONTINUE 
7  =  1+1 

GO  TO  31 

WE  HAVE  DONE  40000  MONTE  CARLO  TRIALS,  JCNT  OF  WHICH  YIELDED 
A  VALUE  OF  THE  INTEGRAND  ABOVE  IIANY.  WE  NOW  CALCULATE  AND 
REPORT  THE  P!CS)  ESTIMATE  AND  IT’S  STANDARD  DEVIATION. 


34  CONTINUE 

i'CSKST°(  JCNT/4Q()C0.)*2.Gl*2.0*!  1 . 0/aQJlT!  2 . 0*3 . 14159)) 
flTDKRa  ( ,)CNT/<'0i?0 ?• .  )*'  1.3-1  JCNT/4U00O. ) ) 

WIDER"  SdllTf  STDER)  /20O . 

STUER=  STDER#2 . 3 1  #2 .  *  (  1 . 0/SCtRT!  2.0*3.14159)) 

WRITE!  6, 01) 

41  FORMAT!  1111) 

WRITE!  6, 42) 

42  FORMAT!  1110, 6X.54HTBE  PROBABILITY  OF  CORRECT  SF.LECTION  IN  A  PROBLEM 

i  wvno 

WIU TEC 0,43)  K 

43  FORMAT!  1II0,6X,5H  K=  ,  12,  12H  POPULATIONS) 

DO  Go  1= I , K 

'WRITE!  0,44)  I ,  XM!  1 )  , SC  I )  , N<  I ) 

44  FORMAT!  UI0.6X,  13H  POPULATION  ,12,1011  HAS  MEAN  ,F10.4, 

HUH,  STD.  DEV.  ,F10.4,19H,  AND  SAMPLE  SIZE  ,14) 

35  CONTINUE 

WRITE' 6, 45)  PCSES-T 

45  FORMAT!  1110, ON,  1 1IIIS  P!CS)  =  ,F6.4,34H  BASED  ON  40000  MONTE  CARLO  T 
1  RIALS) 

WHITE! 0,46)  STREW. 

46  FORMAT!  1R0.6X, 2UILAND  HAS  STANDARD  ERROR  =  ,F6.4) 

WRITE!  6,47) 

47  FORMAT!  till) 

END 
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0601 

0002 

0003 

0004 

0005 

0000 

0007 

000» 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0010 
0017 
00  *« 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
002«'- 
0027 
0020 
0029 
0030 


SUBROUTINE  VSORTC  XM, R, K) 

INTEGER  R, 3 , SAV2 

mmwa  ion  xr«  i  ) ,  R(  i  ) ,  ym(  50) ,  sc  so) 

30  10  1=1,  K 
YHC  I>=XM(  I) 

:*v  i  >  =  i 

10  CONTINUE 

r 

”=K 

1  ;I~n/2 

I  Fell.  EG.  0)GO  TO  3 
J=  t 

2  !  =  .J 

4  1F<  YM-:  I)  .LE.YMC  I+M) )  GO  TO  3 
SAVt=YH(  I) 

i.J  /i  V-s'J'  f  > 

‘>tl< T > •- YMt  IKD 
N\  U=:>(  I+M) 

m  i+vi> =savi 

;v  I-.  K)  -SAV2 
!  =  I-M 

IFl l.GE. 1)  GO  TO  4 

3  .!=«!+ 1 

IF(.J.CT.K-?D  GO  TO  1 
go  TO  2 
3  CONTINUE 
20  11  1=1,  X 
IK  I)  ==  SCI) 

11  continue 

t 
): 


nftfir.nnnnofliionflftof.finnr.nflftnftfinfifif.ftftftftnoft 


FUNCTION  DCDFN(X) 

<t*****«***)(!#*:*#**»**yt*s|i***#***5(«»*****»***C*t:*^**sf:*;f:*******<r*»':5l:*****=l:** 

A  SUBROUTINE  TO  CALCULATE  PROBABILITIES  FOR  STANDARD  NORMAL  RANDOM 
VARIABLES. 

PURPOSE 

GENERATES  T1IE  PROBABILITY  THAT  Z  IS  L ESS  THAN  OR  EQUAL  TO  X, 

WHERE  Z  IS  A  STANDARD  NORMAL  RANDOM  VARIABLE 

USAGE 

“LET  P=DCDFN(  X) * 

DESCRIPTION  OF  PARAMETERS 

X.  .  .  IS  THE  CU*;X>FF  POINT  FOR  WHICH  THE  PROBABILITY  THAT  Z  IS  LESS 
THAN  OR  EQUAL  TO  X  19  DESIRED 

REHAFJCS 

THIS  IS  A  DOUBLE  PRECISION  FUNCTION  AND  CAN  ACCOMMODATE  A  DOUBLE 

hue: i s i on  argument 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

FO.-i  THE  ALGORITHM  'Jo ED  AND  ITS  ACCURACY,  SEE  REFERENCE  2.  FOR 
S  J  ED  AND  C'MPAMSUW  WITH  OTHER  ALGORITHMS  (WITH  REGARD  TO  SPEED 
A SO  ACCURACY* ,  SEE  REFERENCE  I. 

1.  D:J1>EWICZ,K.  J.  AND  HALLEF.T.  :  "TEf> t  XIAND-RANDOM  NUiW.R  GENERATION 
AND  TESTirfC  PACKAGE  <  'NGLUDING  lllCCRAND) ,  “  TECHNICAL  REPORT, 

DEr.Atim-: of  statistics,  the  ohio  state  university,  columbus, 

OHIO.  1  <)'.<.• . 

2.  rm.Ton.n  <:.  ane*  rotcjkiss.r.  :  “con?unsn  evaluation  of  the 

NORMAL  AND  INVERSE  NOii"i\L  DISTRIBUTION  FUNCTIONS,  n 
■JllGlillOITETiiICS .  VOL.  1 1(  1909)  ,  PP.  G 17-022. 

C 

DOUBLE  PRECISION  FUNCTION  DCDFN(X) 
fi  ll  :n;?ION  .\’6>  ,D(0>  ,G(6,  10) 

'  ittAL*0  A . C . D . X. Y . Z  SGN  Y 

•  ATA  v/  .  62:  >1)0 ,  1  .  EU DO ,  2  .  ODD  .  2 . 43D3 , 3 . 5P0 , 4 . 62IV)/ 

DATA  0/.3D).  .92300,  l . 6231>'' , 2 . 223  DO , 2 . 05 DO ,  A .  KiDO/ 

D.vTA  C/f>  .7  ;,)2  i032U  ID-04,-  1 . 270  )7333 93D- 03,0. 7 ‘;ft4N2“797D-04, 


*7 . 27033  i  1 0  ;3D-02 . 4 . 7042  VIR207D-02 ,  -2 .  **•  0347:  iOE/O’D  -f-2 , 

*-  1  ■  022  1  -  29.  ;3  >D-02 ,  -6 .  '  .38  !.V’,0002D-04  .  -4 . 210“-)<K094i)-07, 

*- 1 . 4oooayuiK:4D-o  1,3. 00392  >3gl>6d- 02 , 5; .  74293  isbidd-os. 


Best  Available  Copy 
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*1.  1 05G195&35D-02.5.  12Gn57G121D-04,2.O?42in33<V>I>-O7, 

00011610 

*-  I .  54609 1 297«J»-0 1,-3. 2 1806  1 5 1389-0 1,-6. 53;  13705  1 7OD-02 , 

00011620 

*~a .  «a0403097ttD-  !3 .  -2 . 7637  16253.2.9-04 ,  -7 . 6&“S<>492 1 8D-00, 

0001)030 

.  15630434609-0  1,2.  39790439789-0  1,4. 0236  129 4799-02 . 

00(-«J.  1640 

»r, .  <>93oo.‘?av.  >  *. ,  9 . 375 14 im 4070-03 , ) .  8722020.:;229-g:; , 

00 1  ioro 

a  1.4431 05797  l.)-9> ,  4 . 043033 ’SZ  150-01 ,4.  8922  Tfi  >6599-01 , 

COOl 1660 

*4 . 99  174 167:26>-0  i ,  4 . 9V‘J8489C-*9D-d  1,4. 9999999779!)-© 14 

COOl 1070 

v-  :ii*o .  707  lOvra ;  19  do 

0001  a-m 

OCNY-  1 .00 

0001 1090 

1I'V  Y)8,  1 ,3 

0001  1.700 

1 

TO  T’-  .5D0 

COO 1 1710 

0001 1720 

o 

SGilY-- 1 .  DO 

0001  (709 

V«  _V 

0001. 1 7  40 

3 

<M)  4  1=  1  ,<i 

000  y  ?7;.o 

111  Y- A(  l))5,y,4 

(>0Ol  (760 

4 

O.IMTIWUE 

(-001  1770 

V.- .  UHO 

0001 27150 

I’JO  TO  7 

OOOI  i.7<‘(i 

5 

>.'=  Y-1X  1) 

(■•00;  ir.,/0 

A=C'  1.1) 

0001 10:0 

V*'  <3  .1-2.10 

0001  if  20 

6 

;Vt-C(.  i,.n 

000 1 1  :<'o 

7 

)»ODf’K«  .GDOrUGSYaZ 

(001  KTO 

oooi irro 

liixD 

OOOUG-00 

Best  Ava//ab 
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URK20 

CSECT 

« 

USING 

UR520, 15 

STM 

14,  12, 12(  13) 

ST 

13.SAVE+* 

LR 

11,13 

LA 

13, SAVE 

ST 

13, 8(  11) 

* 

LM 

2,4,«< 1) 

* 

L 

11 ,0(4) 

L 

lO.sF'O' 

L 

*),  -X’ 40000000 

LOOP 

L 

7,0(2) 

M 

4 ,  HLT 

SLDL 

'>*  1 

SHL 

7,  1 

All 

6.7 

* 

C 

n.=F*0’ 

B5L 

UP1 

SLL 

6,  1 

SRL 

6,  1 

4,  =  F’ 1* 

A 

* 

UP1 

ST 

6,0(2) 

LPR 

6 ,6 

L 

u.^x'-fEoeooeo 

STM 

•1,6,  DOUBLE 

LD 

i, DOUBLE 

A!) 

2,=D’O.0’ 

Di) 

2.K0DUL 

STD 

2, 0( 10.3) 

LA 

10, 4(  10) 

BGT 

1 J , LOOP 

L 

13.SWE+4 

LM 

14, 12. 12( 13) 

BCR 

15.  14 

MLT 

DC 

F’ 2027812860’ 

DOUBLE 

DS 

D 

NODUL 

DC 

X’  407FFFFFFFW 

SAVE 

DS 

END 

18F 

ffo jM  UVJiJ' 


;^dii}TOBDC 


COPY  ADDH  OP  CALL  I1TC  SAVE  1*70  PEG  11 
ADDS  OF  SAVE  Ilf  REC13 

STORE  ADDR  OF  SAVE  15  CALL IRC  PROC  SAVE 

REG(2)SADDR(  IX)  REC(  3) * ADDR(  X) 

REGf  4)  =  ADDR(  IfBATCH) 

VALUE  OF  IfBATCH  15  REC11 
IMIT  INDEX 

’  EXFORKMT  FOR  FLOAT!  5C  POI5T  CON  VERS  105 
VALUE  OF  IX  15  REC7 
I.EC6-QU0T,  REC7=REM  >*)D  2**32 
PUTS  OUOT  MOD  2**3I  15  REC6 
PUTS  REM  MOD  2**31  15  REC7 
REG6  =  REII  F»D  2**31-I* 

REM  MOD  2**31  ♦  ODOT  HOD  2**31 
CHECK  FOR  OVERFLOW 
BRANCH  TO  UP!  IF  SUM  IS  ROT  5EG 
D IV IDF.  RIX6  BV  2*^31 
REM  MOD  2*»31  15  REC6 
ADD  ttDOT  FROM  »IV«810R  TO  HEII  TO 
OBTAIN  REM  MW)  2**31- I  15  REG6 
STORK  VALUF.  OF  IX  FOR  NEXT  R5D  bOS 
ABSOL  VALUE  OF  IX  15  REC6 
EXP05K5T  15  REGS  (16**14) 

SHORE  rT.T  PT  50S  15  RECSSJ  fl  6 
LOAD  FLT  PT  REGS  WITH  FLT  PT  50S 
FORMALIZE  F.Y  ADDI5C  O. 

DIVIDE  BY  2**3!-I 

STORE  NORMALIZED  FLT  PT  508  15  X(  I5DEX) 
lilCREKEBT  IRDEX 
BOTTOM  OF  L'M)P 
RESTORE  REGISTERS 

RCTUR5  TO  CALL I 5C  PROGRAM 


=D’0.0’ 

=  F  ’  O  ’ 

-X' 40000000’ 
*F*  1' 

=  X’ 4E0OQOO0  * 


OGCCUOIH 

(  C0«:(.02Ck 
OGC'KW- 
0000' 
(icao-.rW- 

CHXJfn’OAO 

oeoocoro 

OOOOOOoD 


©chxkmo 

00000123 

OOfXX'iGD 

OCOOOIGC 

oooio;  co 

COCO? 160 

OOfOCiTF 

0000017 ( 

0000 '>172 

CC0'’O*7* 

00(  Vvf7<- 

0C(  f  'V  t77 

0C0C'>171\ 

GOOD'  179 

OOtv-ClVf 

OGGOtC-’O 

f  If 

(icoLf:,:v 

ooox:„sf> 

OOCOMK-C 

c-oof*o:*3o 

OOCAW’lfiO 

ooooc-iwi- 
(*cof','>L;y. 
OOOOf LVO 
OOCI'OGJ  0 
000003 1 0 

OOOOOL  4(' 
©COW  l» 
GOGCAoOC 
00000370 


80Sf  Avai,zb!e  Co0l 
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EXAMPLE  AND  COMPARISON 


As  an  example,  when  we  have  k  =  8  populations  with  means  all  0.0  except 
for  one  which  is  2.34,  with  standard  deviations  all  9-0 >  and  sample  sizes 
all  8l,  we  know  (from  the  tables  of  Bechhofer  (1954))  that  the  true 
P(CS)  =  .80.  Our  computer  program  yieldB  the  output 


THE  PROBABILITY  OF  CORRECT  SELECTION  Ilf  A  PROBLEM  WITH 


K=  B  POPULATIONS 


POPULATION 

l 

II  A3 

KEAN 

2.3400, 

8TD. 

DEV. 

9.O0O0, 

AND 

SAMPLE 

SIZE 

81 

POPULATION 

H.’S 

MEAN 

0.0 

» 

STD. 

DEV. 

9.O0O0, 

AND 

SAMPLE 

SIZE 

81 

POPULATION 

y 

HAS 

MEAN 

0.0 

f 

8TD. 

DEV. 

9 . 0000 , 

AND 

SAMPLE 

SIZE 

81 

POPULATION 

■i 

HAS 

MEAN 

0.0 

* 

STD. 

DEV. 

9.0'JOO, 

AND 

SAMPLE 

SIZE 

81 

POPULATION 

(j 

H/TS 

MEAN 

0.0 

« 

STD. 

DEV. 

9 . 0000 , 

AND 

SAMPLE 

SIZE 

81 

POPULATION 

0 

Hlkl 

MEAN 

0.0 

* 

STD. 

DEV. 

9 . 0000 , 

AND 

SAMPLE 

SIZE 

81 

POPULATION 

7 

11A3 

Mi;  AN 

0.0 

■ 

STD. 

DEV. 

9.0000, 

AND 

SAMPLE 

SIZE 

81 

POPULATION 

a 

HAS 

MEAN 

0.0 

* 

STD. 

DEV. 

9.0000, 

AND 

SAMPLE 

SIZE 

81 

IS  PICS)  =  0.8130  BASED  ON  40000  MONTE  CARLO  TRIALS 
AND  HAS  STANDARD  ERROR  =  0.00114 


which  agrees  with  the  theoretical  value,  thus  furnishing  a  check  on  the  computer 


program 
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